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Abstract 

We study the kinetics of chiral phase transitions in quark matter. We dis- 
cuss the phase diagram of this system in both a microscopic framework (using 
the Nambu-Jona-Lasinio model) and a phenomenological framework (using 
a Landau free energy). Then, we study the far-from-equilibrium coarsen- 
ing dynamics subsequent to a quench from the chirally-symmetric phase to 
the massive quark phase. Depending on the nature of the quench, the sys- 
tem evolves via either spinodal decomposition or nucleation and growth. The 
morphology of the ordering system is characterized using the order-parameter 
correlation function, structure factor, domain growth laws, etc. 

Keywords: Chiral symmetry breaking, Ginzburg-Landau expansion, 
TDGL equation, domain growth, quenching, dynamical scaling. 

1. Introduction 

The nature of the quantum chromodynamics (QCD) phase diagram as 
a function of temperature (T) and baryon chemical potential (/z) has been 
studied extensively over the last few years [I!]. For /i = 0, finite-temperature 
calculations have been complemented by first-principle approaches like lat- 
tice QCD simulations jif. However, for fj. ^ 0, the lattice simulations are 
limited to small values of /i [sj. In QCD with two massless quarks, the chiral 
phase transition is expected to be second-order at zero baryon densities. In 
nature, the light quarks are not exactly massless and the sharp second-order 
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transition is replaced by a smooth crossover. This picture is consistent with 
lattice QCD simulations with a transition temperature ~ 140 — 190 MeV 
0] . On the other hand, calculations based on different effective models show 
that the transition becomes first-order at large /x and small T [5|. This means 
that the phase diagram will have a tricritical point (TCP), where the first- 
order chiral transition becomes second-order (for vanishing quark masses) 
or ends (for non- vanishing quark masses). The location of the TCP (/itcp, 
Ttcp) in the phase diagram and its signature has been been under intense 
investigation, both theoretically and experimentally [gI, 0| • 

Heavy-ion collision experiments at high energies produce hot and dense 
strongly-interacting matter, and provide the opportunity to explore the phase 
diagram of QCD. The high-T and low-/i region has been explored by recent 
experiments in the Relativistic Heavy Ion Collider (RHIC). This region will 
also be studied by planned experiments in the Large Hadron Collider (LHC). 
Further, future heavy-ion collision experiments like the Beam Energy Scan at 
RHIC, FAIR in GSI, and NICA in Dubna plan to explore the high-baryon- 
density region of the phase diagram, particularly around the TCP [8|. The 
experiments at the RHIC provide clear signals that nuclear matter undergoes 
a phase transition to partonic phases at sufficiently large values of the energy 
density. However, the nature and kinetics of this transition remains an open 
question. We note here that lattice QCD assumes that the system is in 
equilibrium, whereas heavy-ion experiments are essentially nonequilibrium 
processes. Therefore, information about which equilibrium phase has the 
lowest free energy is not sufficient to discuss the properties of the system. One 
also has to understand the kinetic processes which drive the phase transition, 
and the properties of the nonequilibrium structures that the system must go 
through to reach equilibrium. 

In this paper, we study the kinetics of chiral transitions in quark mat- 
ter. We focus on far-from-equilibrium kinetics, subsequent to a quench from 
the disordered phase (with zero quark condensate) to the ordered phase. 
This rapid quench renders the disordered system thermodynamically unsta- 
ble. The evolution to the new equilibrium state is characterized by spatio- 
temporal pattern formation, with the emergence and growth of domains of 
the preferred phases. This nonlinear evolution is usually referred to as phase 
ordering dynamics or coarsening or domain growth j^, [l3]- Previous studies 
of ordering dynamics in quark matter, which we review shortly, have primar- 
ily focused upon early-time kinetics and the growth of initial fluctuations. 
The present paper is complementary to these studies. We investigate the 



2 



universal properties (e.g., growth laws, scaling of correlation and structure 
functions, bubble dynamics, etc.) in the late stages of chiral kinetics. These 
properties are robust functions of the evolution dynamics, and only depend 
upon general features, e.g., scalar vs. vector order parameters, defect struc- 
tures, conservation laws which govern dynamics, relevance of hydrodynamic 
effects, etc. Our results in this paper are obtained using a time- dependent 
Ginzburg- Landau (TDGL) model, which is derived from the Nambu-Jona- 



Lasinio (NJL) model with two light flavors llj. However, we expect that 



our results apply to a much larger class of systems belonging to the same 
dynamical universality class. 

To place our work in the proper context, we provide an overview of stud- 
ies of dynamical properties of quark matter. These focus on either (a) critical 
dynamics, i.e., time-dependent behavior in the vicinity of the critical points, 
or (b) far- from- equilibrium dynamics, which was explained above. In stud- 
ies of critical dynamics, much interest has focused on the TCP. The static 
universality class of the QCD transition (for non-vanishing quark masses) is 
believed to be that of the d = 3 Ising model, but there is debate regard- 
ing the dynamical universality class. For example, the TCP dynamics was 
argued 12| to be in the class of Model C in the Hohenberg-Halperin classi- 



fication scheme [13|. Essentially, the argument was that critical dynamics is 
described by a nonconserved order parameter (the quark condensate {^pip)), 
in conjunction with conserved quantities like the baryon number density. If 
one includes the mode coupling between the quark condensate and the baryon 
density, the appropriate universality class is Model //with a different dynam- 



ical exponent [IJ]. However, it was also argued that reversible couplings can 
play a crucial role in QCD critical dynamics, which may then differ from that 
of Model H 15 1 . In related work, Koide and Maruyama [l6[ have derived a 



linear Langevin equation for the chiral order parameter. This is obtained by 
applying the Mori projection operator technique to the NJL model. These 
authors study the solution of the Langevin equation, and investigate critical 
slowing down in the vicinity of the TCP. 

Let us next turn our attention to studies of far-from-equilibrium dynamics 



in quark matter. Sasaki et al. [17|] have emphasized the importance of spin- 



odal decomposition in understanding the chiral and deconfinement transitions 
in heavy-ion collisions. Their study was based on a mean-field approximation 
to the NJL model, as well as a phenomenological Landau theory. Sasaki et al. 
discussed fluctuations of the baryon number density as possible signatures of 
nonequilibrium transitions. However, they did not study the corresponding 
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evolution dynamics. Scavenius et al. |18j have investigated the possibihty 
of nudeation vs. spinodal decomposition in an effective field theory derived 
from the nonlinear sigma model. Again, these authors have not investigated 
time-dependent properties, which are of primary interest to us in this paper. 

An important study of evolution dynamics in quark-gluon plasma (QGP) 
is due to Fraga and Krein [l^, who modeled the relaxation to equilibrium 
via a phenomenological Langevin equation. (We follow a similar approach 
in this paper.) This Langevin equation can be derived from a microscopic 



field-theoretic model of kinetics of the chiral order parameter [20|, |21|. Fraga 
and Krein studied the early-time dynamics of spinodal decomposition in this 
model both analytically and numerically, and focused upon the effect of dis- 



sipation on the spinodal instability. In recent work, Bessa et al. [22[ studied 
bubble nucleation kinetics in chiral transitions, and the dependence of the 
nucleation rate on various parameters. 

Skokov and Voskresensky 23| have also studied the kinetics of first-order 
phase transitions in nuclear systems and QGP. Starting from the equations of 
non-ideal non-relativistic hydrodynamics (i.e., Navier-Stokes equation, conti- 
nuity and transport equations), they derived TDGL equations for the coupled 
order parameters. These TDGL equations were studied numerically and ana- 
lytically in the vicinity of the critical point. Skokov- Voskresensky focus upon 
the evolution of density fluctuations in the metastable and unstable regions 
of the phase diagram, and the growth kinetics of seeds. They clarify the role 
of viscosity in the ordering kinetics. Finally, we mention the recent work 



of Randrup [2J], who has studied the fluid dynamics of relativistic nuclear 
collisions. The corresponding evolution equations reflect the conservation of 
baryon charge, momentum and energy. Randrup studied the amplification 
of spinodal fluctuations and the evolution of the real-space correlation func- 
tion and the momentum- space structure factor. (We will study the scaling 
of these quantities in Sec. [3] of this paper.) Randrup's work mostly focused 
upon the evolution in the linearized regime, where there is an exponential 
growth of initial fluctuations. 



This paper is complementary to Refs. [19|, [22, l23|, l2J], and investigates 
the late stages of phase-separation kinetics in quark matter. The system is 
described by nonlinear evolution equations in this regime: the exponential 
growth of initial fluctuations is saturated by the nonlinearity. We study the 
coarsening dynamics from disordered initial conditions, and the scaling prop- 
erties of emergent morphologies. We consider an initially disordered system 
which is quenched to the symmetry-broken phase through either the second- 
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order line (relevant for high T and small /i) or the first-order line (relevant 
for small T and large /i) in the NJL phase diagram. We study domain growth 
for both types of quenching, and highlight quantitative features of the coars- 
ening morphology. We also study the evolution kinetics of single droplets, 
and the dependence of the front velocity on system parameters. 

This paper is organized as follows. In Sec. |21 we discuss the equilibrium 
phase diagram of the two-fiavor NJL model using a variational approach. 
Here, we describe chiral symmetry breaking as a vacuum realignment with 
quark-antiquark condensates. As we shall see, this method also captures 
some extra contributions proportional to l/N^ (where Nc is the number of 
colors), as compared to mean- field theory. In Sec. |21 we will also discuss 
the corresponding Landau description of chiral transitions. In Sec. [21 we 
introduce the TDGL equation which describes the evolution of the chiral 
order parameter, and use it to study the kinetics of chiral transitions. As 
mentioned earlier, we focus on pattern formation in the late-stage dynamics, 
which is characterized by scaling of the evolution morphologies, and the cor- 
responding domain growth laws. Finally, we end this paper with a summary 
and discussion in Sec. HI 

Our investigation has several novel features from the perspective of both 
QCD and domain growth studies. These can be highlighted as follows. First, 
we demonstrate a quantitative mapping between the phase diagrams of the 
NJL model as an effective model of QCD at low energy and the ■i/'^-Landau 
potential. This mapping enables us to identify the relevant time-scales and 
length-scales in chiral transition kinetics. Second, we clarify the quantitative 
features of the coarsening morphology, e.g., correlation functions, growth 
laws, etc., in chiral transitions. These universal features are independent 
of system and model details, and can be measured in experiments on quark- 
gluon plasma. Third, the chiral transition provides a natural context to study 
ordering dynamics in the V^^-potential, which has received little attention. To 
date, most studies of domain growth have focused on the ■i/'^-potential, which 
has a much simpler phase diagram. Finally, as we will discuss elsewhere, 
chiral dynamics also provides a framework to study the effect of inertial terms 
in phase ordering kinetics. Studies of domain growth have almost entirely 
focused on dissipative overdamped dynamics 




5 



2. Phase Diagram for Chiral Transitions 



2.1. An Ansatz for the Ground State 

For the consideration of chiral symmetry breaking, we denote the pertur- 
bative vacuum state with chiral symmetry as |0). We then assume a specific 
vacuum realignment which breaks chiral symmetry because of interactions. 



Let us first note the quark-field operator expansion in momentum space |25 



= ldy'-'[Uo{k)q'i{k) + Voi-k)f^^^^ (1) 

where 

The superscript indicates that and g° are two-component operators which 
annihilate or create quanta, and act upon the chiral vacuum |0). We have 
suppressed here the color and fiavor indices of the quark-field operators. The 
function in the spinors of Eq. ([2]) is obtained as cot0°(A;) = rrii/k 

for free massive fermion fields, i being the fiavor index. For massless fields, 
00 (fc) = 7r/2. 

We now consider vacuum destabilization leading to chiral symmetry break- 
ing [il], described by 

|vac)=WQ|0), (3) 



where 



Uq = exp 



j dk qf{k)\a ■ k)hi{k)qf{-k) - h.c. 



(4) 



Here, hi{k) is a real function of |/c| (= k) which describes vacuum realignment 
for quarks of a given fiavor i. We shall take the condensate function hi{k) to 
be the same {hi = h) for u and d quarks. Clearly, a nontrivial hi{k) will break 
chiral symmetry. A sum over the three colors and three fiavors is understood 
in the exponent of Uq in Eq. (jl]) . 

Finally, to include the effect of temperature and density, we write down 
the state at nonzero temperature and chemical potential where 
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P = 1/T. This is done through a thermal Bogohubov transformation of the 
state using thermo-field dynamics (TFD) [i^, 27|. We then have 



where is 



(5) 
(6) 



Here, 



i3t(/3,/i) 



dk 



In Eq. (171), the ansatz functions 9±{k, /3, /i) will be related to quark and anti- 
quark distributions. The underlined operators are defined in the extended 
Hilbert space associated with thermal doubling in the TFD method. In 
Eq. (JTj), we have suppressed the color and flavor indices on the quarks and 
the functions 9±{k, f3, fi). The ansatz functions h{k), 9±{k,(3,^) will be de- 
termined by minimizing the thermodynamic potential in the next subsection. 

2.2. Minimization of Thermodynamic Potential and Gap Equations 

We next consider the NJL model, which is based on relativistic fermions 
interacting through local current- current couplings. It is assumed that glu- 
onic degrees of freedom can be frozen into point-like effective interactions 
between the quarks. We shall confine ourselves to the two-flavor case only, 
with the Hamiltonian 



(8) 



Here, is the current quark mass. We take this to be the same (m, = m) for 
both u and d quarks. The parameter G denotes the quark-quark interaction 
strength. Further, r is the Pauli matrix acting in flavor space. The quark 
operator ip has two indices i and a, denoting the flavor and color indices, 
respectively. The point interaction produces short-distance singularities and, 
to regulate the integrals, we restrict the phase space to lie inside the sphere 
k < A, the ultraviolet cut-off in the NJL model. 

We next obtain the expectation values of various operators for the varia- 
tional ansatz state in Eq. ([5]). One can calculate these using the fact that the 
state in Eq. arises from successive Bogoliubov transformations. These 
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expressions will then be used to calculate the thermal expectation value of 
the Hamiltonian, and to compute the thermodynamic potential. With 
as defined in Eq. ([1]), we evaluate the expectation values: 



and 



Here, 



A!^(fc,/3,/i) = -[lT(sin2e_ -sin2^+)± (7°cos0i + 
d ■ k sin 0j) (l — sin^ ^_ — sin^ 6'+)] . 



In Eq. (ITT]) , we have introduced the notation = (p^{k) — 2hi{k) in favor 
of the condensate function h{k), which will later prove suitable for variation 
of the thermodynamic potential. 

Using Eqs. fP])- ffTT]) . we can evaluate the expectation value of the NJL 
Hamiltonian in Eq. (|8]) as 



(2vr)= 



-G 



dk 
1 



m cos (pik) + k sin 0(A;) x (l — sin^ 0_ — siv? 
1 



2N, 



■pI 



(12) 



where Np is the number of fiavors. In Eq. 0121) . we have set rrii = m and 
(f)i{k) = (f){k). We have also defined the condensates. The scalar condensate 
is 



Ps 



' (27r)3 



j dk cos (f){k) (l - sin^ 6. - sin^ 9+) , (13) 



and the expectation value of the number density is 



2N^Nf 
(27r)3 



J dk (sin^ e_ - sin^ 9+) . 



(14) 
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The thermodynamic grand potential is then given by 



n = e- fip^- is, (15) 



where s is the entropy density for the quarks. We have the expression [26 



/ dk sin 6. In (sin^ + cos^ 9^ In (cos^ 9^) 

+ sin^ 9+ In (sin^ 9+) + cos^ 9+ In (cos^ 9+) ] . (16) 



Now, if we extremize Q with respect to h{k), or equivalently with respect 
to the function (j){k), we obtain 

M 

cot(/.(A;) = — , (17) 



where M = m — 2gps with g = G{1 + 1/ANc). Substituting this in Eq. f|T3|l . 
we have the mass gap equation for the quarks as 

M = m + 2g'^^^ f dk , ^ (1 - sin^ 9, - sin^ 9+). (18) 

Similarly, minimization of the thermodynamic potential with respect to the 
thermal functions 9±{k) gives 

exp(/^w±) + 1 

where oj± = y/k^ + ± z/ = e{k) ± v. Here, z/ is the interaction-dependent 
chemical potential given as 

G 

= f^- j^Pv (20) 

In Eq. f lT5]) . we substitute the expression for the condensate function from 
Eq. f[T7|) . and the distribution functions from Eq. f lTI?]) . to obtain 



fi(M,/3,/i) = -^-li^y"ci^ {ln[l + exp(-/3u;_)] + ln[l + exp(-/3a;+)]} 

I dkVWTW + gpl-^pl (21) 



(2vr) 



where we have written down the expression for A^^^ = 3 and Np = 2. From 
Eq. ( ETj) . we subtract the potential of the non-condensed state at T = and 
/i = to obtain 



fi(M,/3,/i) 



Q{M, P, fi) -Qoim,P = oo,fi = 0) 

~(2^ / {ln[l + exp(-^u;_)] + ln[l + exp(-/3u;+)]} 
12 



(2vr)= 



j dk (yWTJp - VWTm^^ 



+9Ps - 9Pso - 



G 



pI- 



Here, is the scalar density at zero temperature and quark mass m: 
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(2 



TT 



m 



(22) 



(23) 



2.3. Landau Theory for Chiral Transitions 

In the mean-field approximation and near the chiral transition line, the 
thermodynamic potential obtained above can also be described by Landau 
theory. Let us focus on the case with zero current quark mass. We consider 
the potential in Eq. ( I22l) with m = 0, and terms of order A*""^ being neglected 
(i.e., Nc — )■ oo): 



fi(M,/3,/i) 



12 



(27r)3/3 



dk I In 



+ In 



12 



(2vr)^ 



dk (ywTW -k^ + —. (24) 



To compute this potential numerically, we set the three-momentum ul- 
traviolet cut-off A = 653.3 MeV, and the four-fermion coupling G = 5.0163 x 



30 



10-6 MeV 
and T = is M 
in Fig. [It^a). For fi < fii{T 



With these values, the constituent quark mass at /x = 
312 MeV. The variation of M with /x at T = is shown 
0) ^ 326.321 MeV, the quark masses stay at 
their vacuum values. A first-order transition takes place at fi = fii, and the 
masses of u and d quarks drop from their vacuum values to zero. In Fig. [I](b), 
we show the T-dependence of M at /i = 0. Chiral symmetry is restored for 
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quarks at T ^ 190 MeV. In this case, the transition is second-order: this is 
reflected in the smooth variation of the mass, which is proportional to the 
order parameter (V'^)- 

In Fig. Et^a), we show the phase diagram resulting from Eq. f l24|) for the 
chiral transition in the [fi, T)-plane. The solid line is the critical line, and 
corresponds to the chiral phase transition, which can be either first-order 
or second-order. The first-order line I (at high fi and low T) meets the 
second-order line II (at low fi and high T) in a tricritical point (/itcp, ^tcp) — 
(282.58, 78.0) MeV. A first-order transition is characterized by the existence 
of metastable phases, e.g., supersaturated vapor. The masses corresponding 
to these metastable phases are local minima of the potential, but have higher 
free energy than the stable phase. The limit of metastability is denoted by 
the dot-dashed lines (5*1 and S2) in Fig.[2]^a) - these are referred to as spinodal 
lines. 

Before proceeding, we should stress that the phase diagram in Fig. Wia) 
only considers homogeneous chiral condensates. However, recent calculations 



by Carignano et al. [31|, Sadizkowski and Broniowski [32[, and Nakano and 
Tatsumi jssi show the existence of inhomogeneous chiral-symmetry-breaking 
phases in the NJL model, e.g., domain- wall solitons, chiral density waves, 
chiral spirals. In that case, the first-order line (and the associated spin- 
odal lines) in Fig. [2]^a) may be replaced by second-order transitions between 
inhomogeneous phases. In this paper, we confine ourselves to the kinetics 
of phase transitions between the homogeneous phases in Fig. I21^a). How- 
ever, it is also of great interest to study the ordering dynamics from (say) a 
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Figure 1: (a) Variation of the mass gap M (proportional to the scalar order parameter 
{■ipip)) with quark chemical potential /i at T = 0. (b) Variation of M with T at = 0. 



11 



homogeneous phase to an inhomogeneous phase, or between different inho- 
mogeneous phases. For example, there have been some studies of ordering 
to a lamellar (striped) phase in Rayleiqh-Benard convection, described by 



the Swift-Hohenberg equation [SJ, |35|, |36|, |37|. Another important system 



with ordering to inhomogeneous phases is that of phase- separating diblock 



copolymers [38|, |39|, |40 



Close to the phase boundary, the thermodynamic potential (which is even 
in M) may be expanded as a Landau potential in the order parameter M: 

h{M) = n (0) + -M' + -M^ + -M^ + 0(M«) = fiM), (25) 
2 4 6 



correct upto logarithmic factors 17|, |28|] . In the following, we consider the 
expansion of Cl (M) upto the M^-term. This will be sufficient to recover the 
phase diagram in Fig. [2]^a), as we see shortly. The ffist two coefficients [^{0) 




Figure 2: (a) Phase diagram of the Nambu-Jona-Lasinio (NJL) model in the (^,T)-plane 
for zero current quark mass. A Hne of first-order transitions (I) meets a hne of second- 
order transitions (II) at the tricritical point (tcp). We have (/^tcp,7tcp) — (282.58,78) 
MeV. The dot-dashed hnes Si and S'2 denote the spinodals or metastabihty hmits for the 
first-order transitions. The open symbols denote 4 combinations of {fj,, T) with T = 10 
MeV, chosen to represent qualitatively different shapes of the NJL potential. The asterisk 
and cross denote quench parameters for the simulations described in Sec. 13.21 and Sec. 13.31 
respectively, (b) Plot of ^l{M,j3,^) — f2(0,/3,/z) vs. M from Eq. (|24)). For a particular 
(//, r)-value, we denote the free energy by the same open symbol as in (a). The solid 
lines superposed on the potentials denote the Landau potential in Eq. (|25p with a from 
Eq. (Uni), and 5, d being fit parameters (see Table [T]). 
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and a] in Eq. f l^ can be obtained by comparison with Eq. as 



6 



1 3A^ 
2G 



^dk e {In [1 + e-'^(^-^)] + In [1 + e-^^'^+^^j } 



A 







1 



+ 



1 



(26) 



We treat the higher coefficients in Eq. fl25|) (6 and d) as phenomenological 
parameters. These are obtained by fitting Cl (M) in Eq. (125!) to the integral 
expression for Cl in Eq. (12^ . There are two free parameters in the microscopic 
theory (/i and T), so we consider the M^-potential with fitting parameters b 
and d. For stability, we require d > 0. 

In Fig. [2](b), we plot fl (M) — Q (0) vs. M from the integral expression in 
Eq. flM|) . We show plots for 4 values of (/x, T) as marked in Fig. [2t^a). These 
are chosen to represent qualitatively different shapes of the potential. The 
solid lines superposed on the data sets in Fig. [2](b) correspond to the Landau 
potential in Eq. ( I25l) with a from Eq. ( l26l) . and b, d being fit parameters. The 
values of these parameters in dimensionless units are provided in Table [H 



(/i,T) (MeV) 


a/A2 


b 


rfA2 


A = \a\d/\b\^ 


(311.00,10) 


-1.306x10-3 


0.092 


0.439 


0.067 


(321.75,10) 


3.539x10-3 


-0.101 


0.402 


0.140 


(328.00,10) 


6.431x10-3 


-0.111 


0.396 


0.206 


(335.00,10) 


9.736x10-3 


-0.101 


0.265 


0.255 



Table 1: The coefScients (a, b, d) of the Landau potential in Eq. ([25|) for 4 different values 
of /i at T = 10 MeV. These parameters are specified in dimensionless units of A^, A'^ and 
A~^ respectively, where A = 653.3 MeV. The dimensionless quantity A = |a|(i/|6p will be 
useful in our discussion of the dynamics in Sec. [3] 

The order parameter values which extremize the Landau potential are 
given by the gap equation: 

/' (M) = aM + bM^ + dM^ = 0. (27) 
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Figure 3: Phase diagram for the Landau free energy in Eq. ([25]) in the [6/ (dA^), aj (dA'*)]- 
plane. A hne of first-order transitions (I) meets a hne of second-order transitions (II) at the 
tricritical point (top), which is located at the origin. The equation for I is Oc = 3|6p/ (16(i), 
and that for II is Ac = 0. The dashed hues denote the spinodals S\ and S'2, with equations 
= and ag^ = |5p/(4d). The typical forms of the Landau potential in various 
regions are shown in the figure. The open symbols denote the T)-values marked by the 
same symbols in Fig. HJa) . The cross denotes the point where we quench the system for 
< 0. The asterisk in Fig. [2] (a) corresponds to {a/A?,h,dK^) = (-1.591 x 10-2,8.985 x 
10-^ 7.083 X 10-2) or h/{dK^) = 1.269, a/{dK'^) = -0.225. We do not mark this point in 
the figure as it results in a loss of clarity. 



The solutions of Eq. ( l27l) are 

M = 



= Ml = ^ . (28) 

The phase diagram for the Landau potential in [6/((iA^), a/((iA^)]-space is 
shown in Fig. |3l For 6 > 0, the transition is second-order, as for the M^- 
potential. The stationary points are M = (for a > 0) or M = 0, ±M+ (for 
a < 0). For a < 0, the preferred equilibrium state is the one with massive 
quarks. 
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For 6 < 0, the solutions of the gap equation are 



M = 0, a> Ibl^iM), 

M = 0,±M+,±M_, |6|V(4rf) > a > 0, 

M = 0,±M+, a<0. (29) 

As a is reduced from large values, 5 roots appear at a = |6p/(4(i). However, 
this does not correspond to a phase transition. On further reduction of a, a 
first-order transition occurs at Oc = 3|6p/(16(i). The order parameter jumps 
discontinuously from M = to M = ±M+, where M+ = [3|6|/(4rf)]^/2_ 
The tricritical point is located at 6tcp = 0, a^cp = [cf. Fig. I21^a)]. The 
4 combinations of (/x, T)-values marked in Fig. HJ^a) are identified using the 
same symbols in Fig. [31 



3. Kinetics of Chiral Transitions 

3.1. Dynamical Equation 

Let us next study time-dependent problems in the context of the NJL 
or Landau free energy. Consider the dynamical environment of a heavy- 
ion collision. As long as the evolution is slow compared to the typical re- 
equilibration time, the order parameter field will be in local equilibrium. We 
consider a system which is rendered thermodynamically unstable by a rapid 
quench from the massless phase to the massive phase in Figs. [2](a) or [31 In 
the context of Fig. [3l this corresponds to (say) quenching from a > adh) 
to a < ac{b) at a fixed value of b. (Of course, we can consider a variety of 
different quenches.) The unstable massless state (with M ~ 0) evolves via 
the emergence and growth of domains rich in the preferred massive phase 
(with M = ±M+). There has been intense research interest in this far- 
from-equilibrium evolution jo], [lo| . Most problems in this area traditionally 
arise from materials science and metallurgy. However, equally fascinating 
problems are associated with the kinetics of phase transitions in high-energy 



physics or cosmology j4l|, [42|, [43 



In this paper, we focus on domain growth in QCD transitions, modeled 
by the M^-free energy in Eq. f[25|) . The coarsening system is inhomogeneous. 
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so we include a surface-tension term in the Landau free energy: 

K 



VL[M] 



df 
df 



f{M) + -{vM) 
2 4 6 2V J 



(30) 



In Eq. fl30|) . Q [M] is a functional of the spatially-dependent order parameter 
field M(r), and K measures the energy cost of spatial inhomogeneities, i.e., 
surface tension. 

The evolution of the system is described by the time- dependent Ginzhurg- 
Landau (TDGL) equation: 



d_ 
dt 



M (f, t) 



(31) 



which models the over-damped (relaxational) dynamics of M {f, t) to the min- 
imum of fl [M], i.e., the system is damped towards the equilibrium configu- 
ration [l^. In Eq. ( l3T]) . T denotes the inverse damping coefficient. The noise 
term 6{f,t) is taken to be Gaussian and white, and satisfies the fiuctuation- 
dissipation relation [13i ]: 



{e{r,t)) = 0, 
d{r', t')e{r~^', t")) = 2rT6{r' - r~^')6 {f - t") 



(32) 



In Eq. fl32l) . the angular brackets denote an averaging over different noise 
realizations. Replacing the potential from Eq. fl5U]) in Eq. flHT]) . we obtain 

d 

—M (f, t) = -r (aM + bM^ + dM^) + TKV^M + 9 (f, t) . (33) 
dt ^ ' 

We use the natural scales of order parameter, space and time to introduce 
dimensionless variables: 



M 
r 
t 

e 



MoM', Mo = ^/\a\/\bl 
tot', to = {r\a\)-\ 



(34) 
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Dropping primes, we obtain the dimensionless TDGL equation: 



d 

—M (f, t) = -sgn (a) M - sgn (6) - XM^ + V^M + 6 (f, t) , (35) 
at 

where sgn(x) = x/\x\ and A = |a|(i/|6p > 0. The values of A corresponding to 
T = 10 MeV and /x = 311, 321.75, 328, 335 (in MeV) are specified in Table [D 
The dimensionless noise term obeys the fiuctuation-dissipation relation: 

{0 (r,t)) = 0, 
e{r>, t')e{?\ t")J = 2e 6{f' - r")5 {f - t") , 
, T\h\ 



|a|i/2K3/2- 



(36) 



Our results in this paper are presented in dimensionless units of space 
and time. To obtain the corresponding physical units, one has to multi- 
ply by the appropriate dimensional length-scale ^ and time-scale to. For 
this, we need to estimate the strength of the interfacial energy K, and 
the inverse damping coefficient F. The surface tension can be calculated 
as (J = 

v^(|a|3/7l^l) / dz {dM/dzf. For quark matter^ a is poorly known 



and varies from 10-100 MeV/fm at small temperatures - we take cr ^ 50 
MeV/ fml For T = 10 MeV and fi = 321.75 MeV, we then estimate 
^ = V^/l'A^ 2.8 fm. Similarly, we set F ^ ~ '^T/s, where s is a quantity 



of order 1 |19|. |45||. This leads to to = (r|a|)^ — 2.6 fm/s. 

We study the phase-transition kinetics for two different quench possibili- 
ties. The first case corresponds to high T and low baryon density (/i), where 
the quenching is done through the second-order line (II) in Fig. |2]^a) or Fig. |3l 
The corresponding parameter values are (/x, T) = (231.6,85) MeV [marked 
by an asterisk in Fig. EJ^a)]; or {a/A^,b,dA^) = (-1.591 x 10-^,8.985 x 
10-2,7.083 X 10-2) with A = \a\d/\b\^ = 0.14. 

The second case corresponds to low T and high baryon density (/i), where 
the chiral dynamics can probe the metastable region of the phase diagram. 
This can be achieved by shallow quenching through the first-order line (I) in 
Fig. [2]^a) or Fig. [3l i.e., quenching to the region between I and Si. This case is 
studied using parameter values (/x, T) = (321.75, 10) MeV; or (a/A^, b, dA^) = 
(3.53885 X 10-3,-0.1005344,0.4015734) with A = 0.14. These points are 
marked by a cross in the phase diagrams of Figs. [2]^a) and |3l 
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3.2. Quench through Second-order Line: Spinodal Decomposition 

Let us first focus on tlie ordering dynamics for quenclies through the 
second-order hne (6 > 0) in Fig. [3l For 6 > 0, the chiral transition occurs 
when a < 0. The quenched system is spontaneously unstable and evolves via 
spinodal decomposition jsl, [lo[. The relevant TDGL equation is 



M -M^ - + V^M + e (f, t) 



with the dimensionless potential 

/(M) = - 



4 



6 



The free-energy minima for this potential are 

-1 + V1 + 4A 



M = iM. 



± 



2A 



1/2 



(37) 



(38) 



(39) 



We solve Eq. ( 137|) with A = 0.14 numerically using an Euler-discretization 
scheme. We implement this on a c? = 3 lattice of size A^^ (A^ = 256), with 
periodic boundary conditions in all directions. For numerical stability, the 
discretization mesh sizes must obey the condition 



At < 



Ad + aiAx 



2' 



(40) 



where ai = 4 + (1 — vT+4A)/A. This condition is obtained from a linear 
stability analysis of Eq. fl37|) by requiring riumerical stability of fluctuations 



about the stable fixed points in Eq. f l39l) [46|, |47|. For all results shown in this 
paper, we used the mesh sizes Ax = 1.0 and At = 0.1. We have confirmed 
that this spatial mesh size is sufficiently small to resolve the interface region, 
i.e., the boundary between domains with order parameter — M+ and M+. 
Further, we use an isotropic approximation to the Laplacian term V^M: 



V^M (f, t) 



1 



(41) 



which couples each cell to its 6 nearest neighbors (nn) and 12 next-nearest 
neighbors (nnn). Finally, the thermal noise 6{f,t) is mimicked by uniformly- 
distributed random numbers between f— y4„,A„l. We set = 0.5, corre- 



sponding to e = Al{AxyAt/3 = 0.008 in Eq. (JM]). This noise amplitude is 
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adequate to initiate the growth process from a metastable state on reasonable 
time-scales j48[, as we will see shortly. However, the asymptotic behavior of 
domain growth in both the unstable and metastable cases is insensitive to 



the noise term 46, 49 



In Fig. m we show the evolution of Eq. f l37|) with a disordered initial con- 
dition, which consisted of small-amplitude random fluctuations about the 
massless phase M = 0. The system rapidly evolves via spinodal decomposi- 
tion into domains of the massive phase with M ~ +Mj^ and M ~ — M_|_. The 
coarsening is driven by interfacial defects, which are shown in Fig. |H These 
domains have a characteristic length scale L{t), which grows with time. 

The growth process in Fig. H] is analogous to coarsening dynamics in the 
TDGL equation with an M^-potential 0,|io|, i.e., Eq. ^ with A = 0, which 
describes coarsening in a ferromagnet subsequent to a temperature quench 
from T > Tc to T < Tc. Coarsening in the ferromagnet is driven by kinks, 
with the equilibrium profile Ms{z) = tanh(±z/ a/2). We can use the dynamics 
of kinks to obtain a good understanding of this evolution. The domain scale 
obeys the AUen-Cahn (AC) growth law, L(t) ~ t^/^. Typically, the interface 
velocity v ~ dL/dt ~ 1/L, where measures the local curvature of the 
interface. This yields the AC law. The same growth law has also been 
obtained via a closed-time-path formalism of relativistic finite-temperature 



field theory applied to the NJL model [50 



The pattern morphology in Fig. H] is statistically self-similar in time with 
L{t) setting the scale. The morphology is studied experimentally using the 



order-parameter correlation function [9|, Il0l |: 

C{r,t) = ^jdR (M{R,t)M{R + f,t)^ - (^M{R,t)^ (^M{R + f,t) 

(42) 

or its Fourier transform, the structure factor: 

S{k, t) = jdf e^^-"C (f, t) . (43) 

In Eq. V denotes the volume, and the angular brackets denote an aver- 
aging over independent evolutions. As the system is translationally invariant 
and isotropic, C{r,t) and S{k,t) depend only on the vector magnitudes r 
and k. The existence of a characteristic size L{t) results in the dynamical 
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t = 100 




Figure 4: Kinetics of chiral transition after a temperature quench through the second-order 
line (II) in Figs, ^a) or [31 The d—3 snapshots on the left show the interfaces (M = 0) 
at t = 20,50,100 (in units of to). The defects (interfaces) are kinks between domains 
of the massive phase with M ~ +M+ or M ~ — M+. The snapshots were obtained by 
numericahy solving the TDGL Eq. ([57)1 with A — 0.14, as described in the text. The 
frames on the right show a cross-section of the snapshots at z = N/2. 



20 



scaling oi C{r,t) and S{k,t): 

C(f,t) = gir/L), (44) 
S{k,t) = L'^fikL), (45) 

where g{x) and f{p) are scaling functions which are independent of time. 

Let us demonstrate dynamical scaling for the spinodal decomposition 
illustrated in Fig. HI The statistical results presented in this paper correspond 
to the d = 2 case, and are obtained as an average over 10 independent runs 
with 4096^ lattices. In Fig. [S](a), we plot the scaled correlation function 
[C {r,t) vs. r/L] for 4 different times during the evolution. The length scale 
L is obtained as the distance over which the correlation function decays to 
half its maximum value [C(r, t) = 1 at r = 0]. The data sets collapse onto a 
single master curve, confirming the scaling form in Eq. (jSj). The solid line 
in Fig. [5](a) denotes the analytical result due to Ohta et al. (OJK) [sH, who 
studied ordering dynamics in a ferromagnet. The magnet is also described 
by a scalar order parameter, i.e., magnetization. The OJK function is 

C(f,t) = ^sin-i (^e"'-'/^') . (46) 

(The corresponding result for the case with vector order parameter has been 
obtained by Bray and Puri [52|].) Our correlation- function data is in excellent 
agreement with the OJK function, showing that chiral spinodal decomposi- 
tion is analogous to domain growth in a ferromagnet. 

In Fig. [5]|^b), we plot the scaled structure factor [L~'^S{k, t) vs. kL] for the 
same times as in Fig. ^a). Again, the data sets collapse neatly onto a single 
master curve, confirming the scaling form in Eq. (1451) . The scaling function 
is in excellent agreement with the corresponding OJK function. Notice that 



the tail of the structure factor shows the Porod law 53|, S{k,t) ~ {^-(d+i) 
for /c — )■ oo, which results from scattering off sharp interfaces. For an n- 
component vector order parameter, the system shows the generalized Porod 
law |52|, S{k,t) ~ A;-('^+") for k — )• oo. The scalar order parameter studied 
here corresponds to n = 1 and the relevant defects are interfaces or kinks. 
Higher-order defects arise for vector fields, e.g., vortices or vortex strings 
(ra = 2), monopoles (n = 3), etc. 

In Fig. [6l we plot L{t) vs. t on a log-log scale. The length scale data is 
consistent with the AC growth law, L{t) ~ t^^"^. As stated earlier, this law 
has also been obtained from finite-temperature field theory applied to the 



NJL model 50 . 
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Figure 5: (a) Dynamical scaling of the correlation function [C{r,t) vs. r/L] for chiral 
spinodal decomposition at four different times. The length scale L{t) is defined as the 
distance over which C (r, t) decays to half its maximum value. The different data sets 
collapse onto a single master curve. The statistical data shown in this figure is obtained 
on 4096^ lattices as an average over 10 independent runs, and the correlation function 
is spherically averaged. The solid line denotes the OJK function in Eq. gSl) :5H. (b) 
Dynamical scaling of the structure factor [L~^S{k, t) vs. kL] for the same times as in (a). 
The large-fc region (tail) of the structure factor obeys the Porod law [s^, S{k,t) ^ k^^ 
for k ^ 00, which results from scattering off kink defects. 
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3.3. Shallow Quench through First-order Line: Nucleation and Growth 

Let us next consider shallow quenches from the massless state through 
the first-order line (6 < 0) in Fig. [31 The chirally-symmetric phase is now 
metastable, and evolves to the stable massive phase via the nucleation and 
growth of droplets. From Fig. |3l we see that the first-order chiral transition 
occurs for a < Oc = 3|6p/(16(i). (In terms of dimensionless variables, the 
transition occurs for A < Ac = 3/16.) We consider quenches from a > 
(with M = 0) to a < a^. If we quench to a < 0, the free energy has a 
double-well structure, as shown in Fig. [31 Again, the ordering dynamics is 
analogous to that for the ferromagnet. We have confirmed (results not shown 
here) that the domain growth scenario is similar to that shown in Figs. [IB 

Subsequently, we focus only on quenches to < a < or < A < Ac- 
The appropriate dimensionless TDGL equation is 

— M (f, t) = -M + M^- XM^ + V^M + e (f, t) , (47) 
with the potential 

/ (M) = ^M' - ^M' + ^M^ (48) 

The free-energy extrema are located at 

M = 0, ±M+, ±M_, (49) 

where 



2A 



The extrema at M = 0, ±M_|_ are local minima with f{±M^) < /(O) = for 
A < Ac. 

3.3.1. Bubble Growth and Static Kinks 

Before we study the ordering dynamics of Eq. fH7j) , it is useful to under- 
stand the traveling-wave and static solutions. After all, the phase transition 
is driven by the dynamics of kinks and anti-kinks (eventually 1-d in nature). 
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For the case with 6 > and a < discussed in Sec. I3.2[ the kinks have 
tanh-profiles with small corrections due to the M^-term in the potential. 
We consider the deterministic version of Eq. f H7|) in d=l: 



-M{z,t) = -M + M3-AM5 + — -. 
at oz^ 



(51) 



We focus on traveling-wave solutions of this equation, M {z, t) = M {z ~ vt) = 
M [t]) with velocity f > 0. This reduces Eq. flSTl) to the ordinary differential 
equation: 



+ v- 



dM 



drf ' drj 

Equation (l52l) is equivalent to a 2-d dynamical system 

dM 



0. 



(52) 



dr] 



drj 



vy. 



(53) 



To obtain the kink solutions of this system, we undertake a phase-plane anal- 
ysis. The phase portrait will enable us to identify kink solutions of Eq. flHT]) . 

The relevant fixed points (FPs) are {M,y) = (0,0), (±M_,0), (±M+,0). 
We consider small fluctuations about these FPs: 



M 

y 



Mo + 
+ y. 



(Mo = 0, ±M„, ±M+), 



(54) 



We can linearize Eq. (15^ about these FPs to obtain 



d0 
drj 
dy 
drj 



(1 - SM^ + 5XM^) 



vy 



vy. 



(55) 



The eigenvalues (Ae) which determine the growth or decay of these small 
fluctuations are determined from 



-A, 



(56) 
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or 



-V ± Vv^ + 4a 
^e± = ^ • (57) 

We can combine this information to obtain the phase portrait of the 
system in Eq. (!53|) . In Fig. [TJ we show phase portraits for A = 0.14 (< Ac — 
0.1875), which is the parameter value for most simulations presented in this 
subsection. Figure [Tl^a) corresponds to the case with v = (static solution). 
The saddle connections from — M+ — )• +M_|_ and +M_|_ — )■ — M_|_ correspond 
to static kink solutions. These can be obtained by integrating 



s 



dz 



±V2 



I (Mf -Ml)-\ {Mt - M^) + ^ {Mt - Ml) 



1/2 



(58) 



The static kink profiles for several values of A < Ac are shown in Fig. [Hi 

In Fig. [7](b), we show the phase portrait for v = Vs-, where Vs corresponds 
to the appearance of saddle connections from — M_|. — )■ and +M4. — )■ 0. 
These correspond to kinks traveling with velocity Vs > 0, as shown in 
Fig. nih)- So far, our analysis has been done for the case with f > 0, but it 
is straightforward to extend it to the case with f < 0. In the latter case, the 
portrait in Fig. El^b) is inverted, and the saddle connections (kinks) are from 
-)> -M+ and +M+. 

In Fig. [9l we show the growth of a bubble (droplet) of the massive phase 
(M = +M_|_) in the background of the metastable phase (M = 0). These 
snapshots are obtained by solving Eq. fH7|) with A = 0.14 and 6 = 0. We 
start with an initial configuration of a = 2 bubble of radius i?o > Rc such 
that 

M(r) = M+, r < i?o, 

M(r) = 0, r > Rq. (59) 

Here, Rc{\) is the critical size of the droplet, which diverges as A — )■ A~. In 
Fig. [TUT a). we plot the radius of the droplet [R(t) — Rq] vs. t. These curves 
are linear, showing that the bubble of the massive phase grows at a constant 
velocity vb- In Fig. [TOT b). we plot vb vs. A. The quantity (A — Ac) measures 
the degree of undercooling. Our numerical data is in good agreement with 
Vs, which is obtained from the phase-plane analysis [cf. Fig. d^b)]. 

It is also useful to study the limit A = Ac, where there are three coexisting 
solutions of the free energy in Eq. f H5]) . viz., M = 0, ±M+. The phase 
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Figure 7: Phase portraits of the dynamical system in Eq. ((53|) with A 0.14. (a) Case 
with t; = 0. The saddle connections from — +M+ and +M+ — > — M+ correspond 
to static kink solutions, (b) Case with v = Vs = 0.503, corresponding to the appearance of 
saddle connections from — M+ — )• and +M+ — )■ 0. These correspond to kinks traveling 
with velocity Vs > 0. 
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Figure 9: (a) Growth of a bubble or droplet of the massive phase (M = +M+) in a 
background of the metastable massless phase (M = 0). We show the boundary of the 
droplet at three different times. The innermost circle corresponds to the droplet at time 
t = 60. (b) Variation of order parameter along the horizontal cross-section marked in (a). 
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portrait for v = (static kinks) is shown in Fig. [TTJ Now, there are saddle 
connections from ±M+ — )■ and — ±M+. The corresponding kink profiles 
are obtained as solutions of 



±V2(^lMl-\Mf + ^M!y\ (60) 



dz 

supplemented with appropriate boundary conditions. 




Figure 10: (a) Linear growth of a single bubble with time [R{t) — Rq vs. t] for 4 different 
values of A. (b) Plot of the bubble growth velocity vb vs. A. The circles denote our 
numerical data, while the solid line is obtained from the phase-plane analysis (see Fig. [7]). 



3.3.2. Chiral Transition Kinetics 

Next, we consider the evolution from a disordered initial condition. We 
implemented an Euler-discretized version of the TDGL Eq. (147|) on a c? = 2 
lattice of size A^^. In this case, the mesh sizes must obey the numerical 
stabihty condition: 

At < -T^^^^, (61) 



where 02 = — 4 + (1 + \/l — 4A)/A. We used the mesh sizes Aa; = 1.0 and 
At = 0.1. The other numerical details are the same as in Sec. 13.21 

Recall the phase diagram in Fig. |3l We now focus on the region 6 < 0, 
and consider quenches from high values of a (massless phase) to < a < a, 
or < A < Ac- The massless phase is a metastable state of the potential. The 
chiral transition proceeds via the nucleation and growth of droplets of the 



c 
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Figure 11: Phase portrait of the system in Eq. ((53|) with A = Ac = 0.1875 and v = 0. Now, 
there are saddle connections (or static kinks) from ±A/+ — and — )■ ±Af+. 
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Figure 12: Kinetics of chiral transition after a temperature quench through the first-order 
hne (I) in Figs. [D^a) or [31 The d = 2 snapshots on the left show regions with M = at 
t = 200, 400, 4000 (in units of to)- They were obtained by numerically solving Eq. (|47| 
with A = 0.14. The frames on the right show the variation of the order parameter along 
a diagonal cross-section {y — x) of the snapshots. 
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preferred phase (M = ±M^). This nucleation results from large fluctuations 
in the initial condition which seed bubbles, or thermal fluctuations during the 
evolution. This should be contrasted with the evolution in Fig. HI where the 
massless phase is spontaneously unstable and the system evolves via spinodal 
decomposition. In Fig. [121 we show the nucleation and growth process which 
characterizes evolution. At early times {t = 200), the system is covered 
with the massless phase, with only small bubbles of the massive phase. The 
bubbles grow with time (see Fig. |9]) and coalesce into domains (t = 400). 
The coarsening of these domains is analogous to that in Fig. H] - in the late 
stages of the transition, there is no memory of the nucleation which enabled 
growth in the early stages. 




r/L kL 

Figure 13: (a) Scaling plot of the correlation function [C{r,t) vs. r/L] for the evolution 
shown in Fig. [121 The deviation of data sets from a master curve at early times reflect the 
morphological differences between the "nucleation and growth" and "domain coarsening" 
regimes. At later times, the data sets collapse onto a master curve. The solid line denotes 
the OJK result in Eq. pS)) . (b) Scaling plot of the structure factor [L~'^S{k, t) vs. kL] for 
the same times as in (a). At late times, the tail of the structure factor shows the Porod 
law, S{k,t) ~ k~^ for k oo. 

In Fig. [TST a). we plot the scaled correlation function [C(r, t) vs. r/L] 
for the evolution in Fig. [121 (Our statistical data is obtained as an average 
over 10 independent runs with 4096^ lattices.) The morphological differences 
between the "nucleation and growth" [t = 200, 400) and "domain coarsening" 
{t = 2000,4000) regimes is reflected in the crossover of the scaling function. 
At later times, we recover djTiamical scaling and the master function is in 
excellent agreement with the OJK function. Thus, the late-stage morphology 
in Fig. [121 is analogous to that for spinodal decomposition. In Fig. [TST b). we 
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plot the scaled structure factor [L^'^S{k^ t) vs. kL] at the same times as in 
Fig. [rST a). As expected, the structure factor also violates dynamical scaling 
in the crossover regime. 

In Fig. ??(a), we plot the domain size L{t) vs. t for A = 0.12,0.14,0.15 
on a log-log scale. In contrast with Fig. [6l there is almost no growth in the 
early stages when droplets are being nucleated. The growth process begins 
once nucleation is over, with the onset being faster for lower A (or higher 
degree of undercooling). In Fig. ??(b), we plot L{t) — Lq vs. t — to on a 
log-log scale. Here, Lq (~ 10) is the initial length scale and to is the onset 
time for the different values of A in Fig. ??(a). We see that the asymptotic 
regime is again described by the AC growth law, L(t) ~ (t — to)^^^- 

Let us note here that, converting to physical time scales (i.e., multiplying 
the dimensionless time t by to), the time to reach equilibrium seems to be 
very large compared to the typical life-time r of the fire-ball produced in a 
heavy-ion collision (r ~ 10 fm for the RHIC). Let us recall that to — 2.6 
fm and is proportional to F^^. Larger dissipation (larger F"^) will make the 
equilibration time much larger. Thus, in a heavy- ion collision experiment, 
the system may linger in the QGP phase longer during the fire-ball expan- 
sion, even when the temperature has already decreased below the critical 
temperature. Thus, the critical temperature calculated within equilibrium 
thermodynamic models can be higher than the value that shows up in the 
growth of fluctuations in experiments. However, such conclusions depend 
upon our estimate of F. In principle, F can be calculated within the model 
as has been attempted in Ref. jl6| using the Mori projection operator method. 
However, such a calculation is beyond the scope of the present work and we 



have used the estimate of F in Refs. 19, 45 



4. Summary and Discussion 

Let us conclude this paper with a summary and discussion of the results 
presented here. We have studied the kinetics of chiral phase transitions in 
quark matter. At the microscopic level, these transitions are described by 
the Nambu-Jona-Lasinio (NJL) model. In the NJL model with zero current 
quark mass, there can be either first-order (I) or second-order (II) transitions 
between a massless quark phase and a massive quark phase. The lines I and 
II meet in a tricritical point. At the coarse-grained level, chiral transitions 
can be modeled by a Landau potential with an M^-functional. We have 
shown that there is a quantitative agreement between the NJL free energy 
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as a function of (T, /i) and the Landau potential with an appropriate choice 
of parameters. Near the phase boundary, we can identify the coefficients of 
the Landau energy for different values of {T,fi). However, we often consider 
parameter values far from the phase boundary, and it is more appropriate to 
interpret the Landau coefficients as phenomenological quantities. 

We studied the kinetics of chiral transitions from the massless (disor- 
dered) phase to the massive (ordered) phase, resulting from a sudden quench 
in parameters. We model the kinetics using the time- dependent Ginzhurg- 
Landau (TDGL) equation, which describes the overdamped relaxation of the 
order parameter field (scalar condensate density) to the minimum of the 
corresponding Ginzburg-Landau (GL) free-energy functional. We consider 
quenches through both the first-order and second-order lines in the phase 
diagram. There have been some earlier studies of the TDGL equation in 
this context, as discussed in the introductory Sec. [H However, these have 
primarily discussed the growth of initial fiuctuations in the framework of a 
linearized theory. On the other hand, this paper focuses on the late stages 
of pattern formation where nonlinearities in the TDGL equation play an 
important role. 

For quenches through II, the chirally-symmetric phase is spontaneously 
unstable and evolves into the broken-symmetry phase via spinodal decomposi- 
tion. The evolution morphologies show self-similarity and dynamical scaling, 
and can be quantitatively characterized by the order-parameter correlation 
function or its Fourier transform, the structure factor. The domains of the 
massive phase grow as L{t) ~ t^^'^. 

For deep quenches through I, the above scenario applies again. However, 
for shallow quenches, the chirally-symmetric phase is metastable. Then, the 
system evolves via the nucleation and growth of bubbles or droplets of the 
preferred massive phase. In this case, the early-stage dynamics is dominated 
by the appearance of bubbles. The growth and merger of these bubbles re- 
sults in late-stage domain growth which is morphologically similar to that 
for spinodal decomposition. The correlation function, structure factor and 
growth law show a crossover, having different functional forms in the nucle- 
ation and coarsening regimes. 

Before concluding, it is important to discuss the relevance of these results 
for QCD phenomenology and experiments. In the context of heavy-ion colli- 
sions, we make the following observations. Within the uncertainities regard- 
ing values of dimensional quantities for quark matter (e.g., surface tension, 
dissipation), it is not clear whether the system equilibrates completely within 
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the life-time of the firebalL If the system is nearly-equihbrated, the features 
of the coarsening morphology will be similar for quenches through both first- 
and second-order lines in the phase diagram. However, if the equilibration 
time-scale is much larger than the fireball life-time, the morphology is very 
different for quenches through the first-order line, with the system evolving 
through nucleation of bubbles. Consequences of such a first order transition 
have potential relevance since they imply the existence of a critical end point 
(CEP) in the QCD phase diagram. As a matter of fact, experimental studies 
of such signatures may be more convenient than directly searching for the 
CEP via critical fluctuations. The latter approach has not provided conclu- 
sive evidence of the existence of a CEP, presumably due to the smallness of 
the critical region. 

We also stress that relating our results to heavy-ion collision experiments 
requires information about the source size apart from its life-time. In this 
context, two-particle momentum correlations (i.e., the Hanbury-Brown-Twiss 
or HBT effect in heavy- ion collisions) could be relevant. In such two particle 
correlations, the inverse width of the correlation function in the "out" di- 
rection measures the life-time of the source, whereas the same in the "side" 



direction measures the transverse size of the source [5J|. Thus, the results 
presented here regarding domain growth and nucleation of bubbles can, in 
principle, be investigated through two-particle HBT correlations. This is 
similar to the study of the HBT effect for particle production from inhomo- 
geneous clusters of QGP fluid Finally, the spatial inhomogeneities due 
to the order-parameter evolution could also have measurable effects on the 
spatial distribution and integrated inclusive abundances of various hadrons 
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